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Abstract 

We study a solvable model of two types hard core Bose particles. A complete 
analysis is given of its equilibrium states including the proof of existence of 
Bose-Einstein condensation. The plasmon frequencies and the quantum normal 
modes corresponding to these frequencies are rigorously constructed. In particular 
we show a two- fold degeneracy of these frequencies. We show that all this results 
from spontaneous gauge symmetry breakdown. 
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1 Introduction 



The study of lattice Bose gas models was proposed |T|,0] in order to understand the effect 
of interaction on Bose-Einstein Condensation (BEC). It was argumented that, if the in- 
teraction has a hard core, i.e. each site can be occupied by at most one particle, then the 
model can be mapped onto a quantum spin-^ XY-model on a Z'^-lattice. 

The exact solution, exhibiting condensation, for the hard core Bose gas on a complete graph 
has been obtained P,^. In spin language this work corresponds to the mean field spin-i 
XY-model. Other exact results on the full model are an upper bound on the condensate 
0] and an existence proof of condensation for the system in a half filled Z'^-lattice for 
dimensions (i > 3 [0. In this context we mention also the relevant result |^ about the 
occurrence of long-range order in the XY-model with spins and dimensions greater than 
one. 

In this paper we turn our attention to a system of two types of interacting hard core Bose 
particles on a complete graph and interacting through a repulsive interaction if the parti- 
cles are of different type. (Compare this with the Hubbard model for electron systems.) 
The model is exactly soluble. We study in detail the equilibrium states in the case of a half 
filled lattice and equal density for the two particle types. If the system is weakly interact- 
ing we obtain an explicit and complete analysis of the equilibrium states. We analyse the 
BEC phenomenon and in particular we are interested in the appearance of spontaneous 
symmetry breaking (SSB) linked with it. The interesting aspect of this model is that the 
broken symmetry is two-dimensional. While, as is well known for a single type Bose sys- 
tem, the broken symmetry is the one-dimensional gauge symmetry. We compute rigorously 
the plasmon frequencies related to these modes, and find the quantum normal modes of 
these plasmon motions. We find out that there are two different frequencies, but that both 
frequencies are two-fold degenerate. The four quantum degrees of freedom have normal 
modes which we build up as the quantum fluctuation operators of the components of the 
infinitesimal generators of the spontaneously broken gauge symmetries, with adjoints the 
quantum fluctuation operators of corresponding order parameter operators. The degener- 
acy of the plasmon spectrum is the new phenomenon in this model study of Bose-Einstein 
condensation. Another revealing fact is that the adjoint fluctuation operators are the total 
and the relative current fluctuations between the two types of bosons. 



2 The Model and its Equilibrium States 
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2.1 The Model 



The model is defined on the lattice Z, by the local Hamiltonians with A a finite interval 
in Z: 

' ' ae{l,2} i^ieA oe{l,2} ieA ieA 

The first term is the hopping contribution, the particles jump fireely from one site to another 
with a hopping rate t/|A|, |A| is the number of points in A; the index a specifies the type of 
particles; /i is the chemical potential; U is taken to be positive and is the coupling constant 
for the interaction between the types of particles. If two particles of different type are on 
the same lattice point the potential energy is raised with an amount U . Clearly the a^, 
are the usual Pauli matrices and are here the creation an annihilation operators for the 
Bose particles with hard core, satisfying: 

(^i^)^ = 0; (^i)* = f^i;; 

C^ia^^ia = 2 + l) • 

We look first for the equilibrium states of this model in the thermodynamic limit. The 
Hamiltonians (|l|) are permutation invariant with respect to the lattice index. This means P] 
that the equilibrium states at inverse temperatures (3 are convex combinations of product 
states on the infinite tensor product algebra of observables M. = <S>iez (^2 ® ^2)4 of the 
system; (M2 ® ^^2)4 represents the tensor product of the 2x2 complex matrices M2 with 
itself. Each set M2 is referring to a different type of hard core bosons, both at site i G Z. 
In particular, if one looks only at the product states which are invariant for the interchange 
of the two types of particles, then the equilibrium states are given by the states ujpx^ . ) of 
the form P]: 

n 

ujfsx{AiA2 ...An) = Yltipf^xA; for any A, G (M2 ® Ms)^ ,i = l,...n 

i=l 

where ppx is a density matrix in M2 ® M2, satisfying the self-consistency equation: 

- trexp[-/3/..]' 

where 

hx = tX [erf + 0"^) + tX ((Tf + (T2) — /U (o"^(T]~ + 0-2^2) + f^^i'^r^2''^2' (3) 

and 

A = tl Pi3\(Ti = tTpi3X(T2- (4) 
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In fact equation (0) for pp\ is equivalent to the equation (||) for A. Remark that (H) also 
results from the assumption of symmetry in the particle types. Looking for equilibrium 
states breaking this symmetry would violate (^). 

Moreover, we will limit here our study to the case of a half filled lattice. Mathematically 
this means that we require: 

ujpx{crl^) = for alH G Z, a = 1, 2. 

This fixes the chemical potential fi = U/2. In order to simplify the notations, we adapt 
our Hamiltonian by adding a constant term and by scaling the constant U to U/A. In this 
case one gets for a more simple expression: 

hx = t\ {at + at) + t\ ((Tf + da") + Ua^a'^. (5) 

On a finite subset A C Z we define Fourier-transformed operators: 

L |A| . 



A; G A* the dual of A, 



the annihilation operator of the a-particle with quasi-momentum k. 
Then, 

2 

"^k = 0,*ka'^ka 
a=l 

is the number operator of particles with quasi-momentum k. 
Compute in a state Uisx (Q): 



Po = lim uji3x(^) = 2ujpx{(^^)i^px{(^ ) = 2|A|^; (6) 

|A|^oo V|A|/ 



i.e. the 0-mode has a density po equal to 2|Ap. A solution A 7^ of (^) yields a macroscopic 
occupation of the 0-mode and the density of the condensate po equals 2|Ap. Therefore, if 
one finds solutions of (0) or (H) such that A 7^ 0, one proved Bose-Einstein condensation. 



2.2 Equilibrium States 



Now we proceed to look for the solution of and its properties. It is clear that based on 
general results we reduced our problem to a one-site problem, namely looking for the 
density matrix state tr ppx ■ with p^A G M2 (8> M2. 

In general M2 ® M2 = M4 and we use the following isomorphism: 

/ aiifoii auhi 



an 

^21 



ai2 



hi 

hi 



h2 
h2 
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021611 
011621 



O22611 
O'uhi 



aiib 
6^21612 

0'uh2 



012612 \ 
0-22612 
O12622 



\ O21621 a-22621 021622 a-22622 / 
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This yields the following matrix representation for hx: 



( ^ 


tx 


tx 


\ 


tx 


-u 





tx 


t\ 





-u 


tx 


\ 


tx 


tx 


u ) 



(7) 



Then a straightforward diagonalisation of h\ yields the following spectrum: 

eo = -f/, ei = f/, 62 = ?7, es = -r/, 



where t] = ^/U^ + 4|tAp. 

The corresponding eigenvectors are: 



I0o) 



102) 



1 

71 



/ \ 



V / 



V2\tX\ 



/ tx \ 

v-u 



/ tx \ 




\-txj 
/ tx \ 

v+u 



(9) 



Vviv + u) 



v+u 



Using (H) and (|) one gets for equation the following explicit and nonlinear equation 
for A: 



A 



1 + 



t 



sinh(/??7) 



r] cosh(/5f/) + cosh(/??7) 



0. 



(10) 



Remark that A = is always a solution. It corresponds to a state without condensation 
(see m). 



Looking for solutions A 7^ amounts to look for a solution t] satisfying: 
where 



cosh(/5[/) + cosh(/?r7) 



Since the condensate density po = 2|Ap (cfr. (|)) cannot exceed the total particle density 
p (which is equal to one in the half filled lattice), \/lP + 2^ is an upper bound on 77 = 
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^/lPT4\tW■ The minimum value of rj is U. This defines /, the range of accessible values 
for 7]-. 



riel= U, + 2^2 . 

Remark that fi3{i]) is unaffected by the sign of U. On the other hand, equation ([Tl| ) has 
no solutions if t > 0. Therefore we take t < and U > 0. 

Discussion of the solutions of j^ldj) 

We are looking for solutions yielding a second order phase transition. 
Remark first that for r] E I: 

^ > (12) 
or] 



dfM 
df3 



< (13) 
> (14) 



i.e. Tj — > fpiff) is a monotonically increasing concave function on /; also (3 —>■ fi3{i]) is 
monotonically increasing for T decreasing. 

The problem here is that, due to the presence of the potential U, these properties are not 
valid outside the interval J, such that the situation is more complicated than for the usual 
mean field case. 

First we consider temperature T > 0, and we distinguish two cases: 
case a): 

fp{U) > U (15) 

yielding 



2U \-t-2U, 
and implying the condition on the interparticle interaction: 

2U < -t (17) 

i.e. the interparticle interaction must be relatively small compared with the hopping term 
(|^); this condition is a meaningful situation for this kind of systems, if the interaction 
becomes to strong, the typical quantum effects like BEG disappear. 
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One derives also 

fp + 2t2 j < + 2t2. (18) 

Together, condition (^) and ( ]T8| ) yield a unique solution of //^(r]) = r] E I, and because 
of (|l^ a second order phase transition with A 7^ as an order parameter. Equation ( |I6D 
learns that there exists a critical jSc determinated by: 

^ i,„f4±|^V (19) 



ksTc 2U \-t-2U 

This defines the critical temperature Tc. 
case h): 

fm < u (20) 

This case is somewhat more complicated, a careful analysis learns that apart from a second 
order phase transition, a first order one can occur, even if (|17D is satisfied. In order to 
avoid non-natural first order transitions, one has to impose: 



<1; fpXU) = U. 

u 



drj 

The analysis of these yields a stronger condition on U than (|T^): 

U < i-t)K, 

where n can be computed numerically as the solution of the equation: 

, ,'1 + 2k\ Ak 
m 



1-2^7 1-2k2' 
K ^ 0.461292. 

Finally we consider the ground state (T = 0) situation. The self-consistency equation (|T0]) 
becomes: 

7] = lim fp{ri) = f^{r])] 



with: 

fooiv) 



ii7]<U 
-t if r/ > U. 



If U > —t, then there are no non-trivial solutions A 7^ 0. On the other hand, if ?7 < —t, 
then there exists always a unique solution A 7^ 0. It is obtained as the solution of: 

-t= v/f/2 + 4t2|A|2 



7 



given by: 



\\\ = -ie-u^Y"<-. 



Since fp{rD is increasing with T decreasing (0), this also imphes an expression for a 
supremum on the condensate density po 

(t^ - U"^) 

Po < < p. (21) 

The condensate density is always strictly smaller than the total particle density, which is 
compatible with former results on this topic This finishes the study of the equilibrium 
properties of the model (|lD. 



3 Spontaneous Symmetry Breaking 

The model H\ (|l|) and its equilibrium states have a discrete symmetry consisting of the 
interchange of the two particle types into each other. It is implemented by the following 
unitary matrix: u = aia^a2(T2+^i^t^2^2^^i'^2^'^i^t) then for any local observable 
Ai2 depending on the types 1 and 2, one gets: 

71-2(^2) = I \{iUl2)i ) A12 I \{{Ui2)i ) = ^21- 

The eigenvectors of h\ satisfy 

M12I0O) = -|0o); M12I01) = M12I02) = 102); M12I03) = l^s)- 

A second discrete symmetry is the particle-hole {p ^ h) symmetry. This symmetry inter- 
changes the creation and annihilation operators {af and cr") and is implemented by the 
unitaries — '^i'^2 i'^^ the x-component Pauli matrix). The invariance of H\ (|l]) is 
expressed by: 

%^k{Hj,) = ( hks). ) ^A ( nK2)^ 

VieA / \ieA 

The eigenf unctions @ of h\ satisfy 

f^?2l0o) = -|0O); ^?2l0l) = -|0l); ^?2l02) = 102); ^^5^2103) = 103). 
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The Hamiltonian H\ (|1|) has also a two-dimensional continuous gauge symmetry group 
given by the actions: 



l<t>.<t>M) = ^'^"(^t^ e [0,27r], k=l,2. (22) 
For a finite volume A, it is implemented by the unitaries: 



1 



where 



f4(0i,</.2) = exp - ( 0i(5i,A + 4>2Q2,A 



jeA 



Clearly, 'y^^^^i^A) = Ua(4,^,4>2)HaUI(^^^^^,^-^ = Ha- Hence, per lattice site one has a com- 
mutative two-dimensional Lie-algebra of this symmetry group, generated by af and 
However, the equilibrium states up\ with A 7^ (see section H) break this symmetry down. 
Indeed we have here for k = 1,2: 



^/3A (70102 (^fc )) = e '*'=^/3A(fT^) ^ UJ/sxi 



In contrast to all other models for which we have been studying the Goldstone phenomenon 
§,[1^, here we have the spontaneous symmetry breakdown of a two-dimensional group. 



In the following, we construct again the canonical normal modes of the Goldstone bosons 
together with their corresponding plasmon frequencies. We take a symmetry breaking 
state ujf^x with A 7^ and we consider fiuctuation mode operators. Before constructing the 
relevant canonical normal coordinates of the Goldstone particles in the next section, we 
introduce here the notion of fluctuation operators. 



For the mathematical details we refer to [jTT],[T2|. Denote by M|" = (M2 M2y"' the 



vectorspace of the Hermitian 4x4 matrices, i.e. these are one-site observables. 

For any A G M^", its fiuctuation operator Fa{A) for the volume A in the state u^x is as 

usual defined by: 

V M^l xeA 



where is a copy of A on site a: G Z. In [|TT| it is proved that the limit: 



^mF^AA) ^ Fx{A) 

exists in the sense of a non-commutative central limit as an operator on a Hilbert space Tix, 
generated by a normalised vector fix, aiid vectors Fx{Ai) . . . Fx{An)flx, with Ai G Ml", n 
is arbitrary and a scalar product: 
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for Ai,Bj G M|", 

fe(Ai) . . . F^{A^)Qx, F^{B,) . . . F^{BJQ^ = 5„^Perm. ({Q^, F^{A,)F^{B,)Qx)) 

V / \ / i.j=l,...n 

(23) 

where the two-point function is given by: 

(q^, F^{Ai)F^{B,)Q^^ = U3^{AB) - u^^iA) up^{B). (24) 

i.e. by the truncated two-point functions of the given state Ufsx. Furthermore, as aheady 
could be understood from (p^), the fluctuation operators Fx{A) , Fx{B) , . . . satisfy the 
canonical commutation relations: 

[Fx{A),Fxm=^pxi[AB])l (25) 

i.e. the fluctuation operators are quantum canonical variables (compare [q,p] = ih) with 
quantisation parameter tui3x{[A, B]) (instead of ih) and the state (^x, ■ ^xj is a generahsed 
free or quasi-free state of a boson field algebra of observables. In section ^ we look for 
the relevant operators A,B,... in order to obtain the normal canonical variables for the 
Goldstone particles, i.e. the particles generated by the spontaneous symmetry breaking. 
We derive the dynamical equations for these modes and we solve them. 



4 Spontaneous Symmetry Breaking and Plasma Os- 
cillations 

We will consider fluctuations of the generators 

Q± = al± 

of the broken symmetry. 

Let Pi, {i = 0, 1, 2, 3) be the spectral projections on the eigenf unctions corresponding 
to the eigenvalues (||) of hx (|3), i.e. hx = J2^=o^iPi- Consider for any local observable 
A G the Hermitian operators: 

E,,{A) = P,AP, + P.APf, i<3- ^,j = 0,1,2,3; (26) 

and the adjoints: 

JE,, {A) = i {P:,AP^ - P,AP,) . (27) 

An immediate calculation leads to: 

[hx, EkiiA)] = -i(e, - ei)JEkiiAy, 

[hx,JEM{A)] = i{ek-ei)EM{A). ^ ' 
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Now we take A = Q^, respectively A = Q^. Due to the particle-hole and type-exchange 
symmetries (section the only non-zero operators of the type (|26| - p7D are: 



EoriQ.), JEo.,{Q-), EuiQ+), JEuiQ+) 2 = 2,3. (29) 

It is interesting to write out these operators in terms of the Pauli matrices af,a-, i = 
1,2. As an example we give the expression for the first pair of these operators 
(£"02 (<5-), JEo2{Q-)), the others are left as an exercise for the reader: 

Eo2iQ-) = + - oloX) + —{(y\o2 - ola^Y 

zr] rj 7] 

JEq2{Q-) = ^7ri(^i (^2 - (^2(^1) + i— (c^i^ -(^2)- i— (o"r - (^2)- 
irj r] rj 

Eq2{Q-) is a component of the infinitesimal generator Q_ and JEq2{Q-) is the adjoint 
order parameter operator to Eq2{Q-)- 

It is clear from (^) that the dynamical equations ( p8D for the operators (pOj) contain only 
the relative energies — ei given by: 

^± = v±u. 

In particular, the energy gap ^+ is entering in the equations: 

[hx, Eo2{Q-)] = i^+JEo2{Q-), [hx, JEM-)] = ~i^+Eo2iQ-), (30) 
[hx, ErsiQ^)] = ~i^+JE^siQ+), [hx, JE,,{Q+)] = ii+E,,{Q^)- (31) 

and the energy gap ^_ in: 

[hx, Eo;{Q-)] = -i^- Ji?o3(Q-), [hx, JEM-)] = iC-Eo3{Q-), (32) 
[hx, EuiQ+)] = i^-JEi2iQ+), [hx, JEM+)] = -i^-Ei2iQ+). (33) 

All these results are on the level of the microscopic observables. In particular one 
discovers in ( p0| - |33| ) that the dynamics leaves invariant the pairs {Eo2{Q-), JEo2{Q-)), 
{E,s{Q+)^ JEi3{Q+)), (EosiQ-), JEo3{Q^)) and {Eu{Q+), JEM+))- On the other hand 
the different pairs do not form canonical pairs of variables. Therefore, we have to consider 
the macroscopic variables, namely the fluctuations of these quantities. We know from 
section ^ that the fluctuation operators are canonical variables, satisfying the canonical 
commutation relations. 



We treat explicitly the first pair appearing in the equation (|30|) , and give the results for 
the other pairs without technical details. 

Plasmon frequency ^4. 
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With the notations of section ^, denote the fluctuations: 

Xo2 = -—Fx {EM-)) ; Po2 = —Fx {JEM-)) , (34) 

no2 no2 



where 

no2 = ^ J ; Z = 2 cosh(/3[/) + 2 cosh(/3r/). 

These satisfy (following (pSD ) the canonical commutation relation: 

[Xo2, P02] = -^ujpx {[Eo2{Q~), JEM-)]) 1 = (35) 



with quantisation parameter: 



4£ 

n+ = — tanh (/?^+/2) . (36) 

7] 



Following (p^ , the calculations for the variances are straightforward: 

(^x, Xl^nx) = (f^A, Po'2^^a) = 2^. (37) 



Now compute the dynamics of these fluctuation operators induced by the microdynamics 
determined by the Hamiltonian hx (^). In differential form the dynamics follows immedi- 



ately from the equation (|30D : 



^ {Fx {EM~))) = Fx {[hx, Eo2{Q-)]) = i^+Fx {JEM-)) ; 
{Fx {JEM-))) = Fx {[hx, JEo2{Q-)]) = ~i^+Fx {Eo2{Q-)) • 



Using the notations (Q) and in integrated form one gets: 

Xo2{t) =Xo2(0) cos(e+t) + Po2(0) sin(e+t); 



Po2{t) = - Xo2(0) sin(e+t) + Po2(0) cos(e+t). 



i.e. the pair of macroscopic variables {Xq2,Po2) is a canonical pair (|35D , evolving in time, 
independently from all other degrees of freedom of the system, as a pair of quantum har- 
monic oscillator variables oscillating with a frequency ^_|_. Remark that if the temperature 
T tends to the critical temperature T^, i.e. when |A| — > 0, then 0. It follows from the 

equations {^7\j that the pair disappears completely i.e. Xq2 = P02 = for T > Tc. 

This finishes the study of the first pair of fluctuation variables. 
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In a completely analogous way one defines the second pair of canonical fluctuation 
tors: 



with 



Xi3 = — F,(Ei3(Q+)); 
ni3 



ni3 



Pn = —FxiJE,3iQ+)) 
ni3 



One gets, comparable to (|35|), (|37|) and (p8|): 

[Xi3,Pi3] =in^ 



^A, -PlS^A 



Xi3(t) =Xi3(0) cos(e+t) + Pi3(0) sin(e+t); 
P,3{t) = - Xi3(0) sin(e+t) + Pi3(0) cos(^+t). 



Plasmon frequency 



Following (32), consider the third pair of canonical fluctuation operators 



where 



Xo3 = —F,{Eo3{Q^)); 
no3 



no3 



P,, = ^F,{JEo3{Q-)) 
no3 



One gets 



with a new quantisation parameter. 



[Xo3,Po3] =i^- 



4e 



h_ = — tanh i(3^^/2); 



variances, 



^A,-^03^Aj — (px^ Po3^\^ — 
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and dynamics 

Xo3(t) =Xo3(0) cos(e-t) + Po3(0) sin(e-t); 
PosW = - Xo3(0) sin(e-t) + Po3(0) cos(e.t). 

Finally, following (|33D, consider the fourth pair of canonical fluctuation operators: 

Xi2 = -—Fx {Eu{Q+)) ; Pi2 = —Fx {JEu{Q+)) , (48) 

ni2 ni2 

where 



ni2 

One computes the commutator: 



Z 



[Xi2,Pi2] = i^-; (49) 



the variances, 



f2A, X^^r];,) = (^l];,, P2^(];,) = 2^; (50) 



and the dynamics 



Xi2(t) =Xi2(0) cos(e^t) + Pi2(0) sin(e-t); 
Pi2(t) = - Xi2(0) sin(e_t) + Pi2(0) cos(e-t). 



(51) 

This completes the construction of the plasmon spectrum together with all its quantum 
normal modes. The next question is: are these four pairs (0,^,^,^) all different ? The 
question is relevant, because the mode variables are all fluctuations and coarse graining 
might make some of them equivalent, i.e. although A,B & Ml"- and A B, Fx{A) might 
be equal to Fx{B). The problem of equivalence of fluctuation operators is considered before 
(see [0[T^). It turns out that they are equivalent, i.e. Fx{A) = Fx{B), if and only if 

(Qx, Fx{A - BfQx) = ujpx {{A-B- upx {A - B)f) = 0. (52) 
If we take for A, B the operators Eij{Q±) and JEij{Q±) with i < j, then clearly 

ujf3x{Eij{Q±)) = ujpx{JEij{Q±)) = 

and we have to check 

upx {{A - Bf) = ujfsx {A^ + B^-AB- BA) = 0; (53) 
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for A and B taken out of different pairs (i?02(Q-), J Eq2{Q -)) , (-E'i3(Q+), '^-^i3(Q+)); 
{Eq2,{Q -) , J EQ'i{Q J)) , {Ei2{Q+), JEi2{Q+)). It is immediately checked that in this case 
always W/ja i^B) = ujp\ {BA) = 0. Hence, (|53|) reduces to 

oo^x {{A - Bf) = LOfsx {A^ + B^) ^ 0; 

which follows from (p7|) , (PD , (^61) and (pOl) . Hence, the four pairs of canonical fluctuation 
operators are all independent, proving the two- fold degeneracy of the plasmon frequency 
and 

Finally, remark that if the temperature T tends to the critical temperature, i.e. when 
A — s> 0, then the pairs (X02, -P02) (0) and (X13, P13) (^) disappear completely. On the 



other hand the pairs (X03, Pqs) dH) and (X12, P12) (|48|) become classical variables (see ( ^i)) 
and (P^D), and become constants of the motions (see ( ^71 ) and (|^)). In fact for T > the 
continuous symmetry group (^2]) is restored. The fluctuation operators X03 and X12 are the 
fluctuations of the generators of this symmetry. The fluctuation operators Pqs and P12 are 
the total particle current and relative particle current fluctuations. For T < Tc the physical 
interpretation of these operators remains the same with the basic difference that they show 
a definite collective or macroscopic quantum character, and a non-trivial oscillatory time 
behaviour. We remind that the results derived for this model are mathematically rigorous. 
In view of the recent interest in interferences between Bose-Einstein condensates with 



Josephson-like coupling [|T^,|14[ an extension of our model including such a coupling raises 
over as the natural thing to do next. We leave it for a future contribution. 

Acknowledgement: The authors thank Tom Michoel for helpful discussions. 
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